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Abstract 

We compute the leading-order collisional energy loss of a heavy fermion propagating in a QED 
plasma with an electron distribution function which is anisotropic in momentum space. We show 
that in the presence of such anisotropies there can be a significant directional dependence of the 
heavy fermion energy loss with the effect being large for highly-relativistic velocities. We also 
repeat the analysis of the isotropic case more carefully and show that the final result depends on 
the intermediate scale used to separate hard and soft contributions to the energy loss. We then 
show that the canonical isotropic result is obtained in the weak-coupling limit. For intermediate- 
coupling we use the residual scale dependence as a measure of our theoretical uncertainty. We also 
discuss complications which could arise due to the presence of unstable soft photonic modes and 
demonstrate that the calculation of the energy loss is safe. 

PACS numbers: 11.15Bt, 04.25.Nx, ll.lOWx, 12.38Mh 
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I. INTRODUCTION 



An understanding of the production, propagation, and hadronization of heavy quarks in rel- 
ativistic heavy ion cohisions is important for predicting a number of experimental observables 
including the heavy-meson spectrum, the single lepton spectrum, and the dilepton spectrum. The 
first experimental results for the inclusive electron spectrum have been reported ^ in addition to 
the first measurements of J/ ip production at RHIC . The measurement of the inclusive electron 
spectrum allows for a determination of heavy quark energy loss since it is primarily due to the 
semi-leptonic decay of charm quarks. The heavy fermion energy loss comes into play since it is 
necessary in order to predict the heavy fermion energy at the decay point. It is therefore important 
to have a thorough theoretical understanding of heavy fermion energy loss for a proper comparison 
with the experimental results. In this paper we will show that in QED there is a modification of 
the leading-order (collisional) heavy fermion energy loss if there is a momentum-space anisotropy 
in the electron distribution function. The motivation for this work is to provide a testing ground 
for the techniques which can be applied to QCD in order to make predictions of the directional 
dependence of the collisional energy loss of a heavy fermion propagating through an anisotropic 
quark-gluon plasma. 

In the last few years a more or less standard picture of the early stages of a relativistic heavy 
ion collision has emerged. In its most simplified form there are three assumptions: (1) that the 
system is boost invariant along the beam direction, (2) that it is homogeneous in the directions 
perpendicular to the beam direction, and (3) that the physics at early times is dominated by gluons 
with momentum at a "saturation" scale Qs which have occupation numbers of order l/a^. The 
first two assumptions are reasonable for describing the central rapidity region in relativistic heavy- 
ion collisions. The third assumption relies on the presence of gluonic "saturation" of the nuclear 
wavefunction at very small values of the Bjorken variable x j^. In this regime one can determine 
the growth of the gluon distribution by requiring that the cross section for deep inelastic scattering 
at fixed does not violate unitarity bounds. As a result the gluon distribution function saturates 
at a scale Qs changing from l/k'j_ —>■ log{Q'^/k'j_)/as. Luckily, despite this saturation, due to the 
factor of 1/as in the second scaling relation the occupation number of small- a; gluonic modes in 
the nuclear wavefunction is still large enough to determine their distribution function analytically 
using classical nonlinear field theory In the weak-coupling limit the assumptions above have 
been used by Baier et al. in an attempt to systematically describe the early stages of quark-gluon 
plasma evolution in a framework called "bottom-up" thermalization 

The resulting picture which emerges from using these assumptions is one in which the initial 
gluonic distribution function is extremely anisotropic in momentum space having the form 

/(p,x) = F(pxM(p.). (1) 

This is, of course, an idealization. In a more realistic scenario the delta function above would 
have a small but finite width which increases as a function of time, but despite this finite width, 
the distribution function would still be extremely anisotropic in momentum space during the early 
stages of the collision. In anisotropic systems it has been shown that the physics of the QED 
and QCD collective modes changes dramatically compared to the isotropic case and instabilities 
are present which can accelerate the thermalization and isotropization of the plasma 0, 0, 0, 
H, 0, 0, 0]. In fact, the paper of Arnold, Lenaghan, and Moore points out that the presence 
of anisotropics in the early stages of QGP evolution requires modification of the "bottom-up" 



2 



thermalization scenario 




In our previous paper [lO| we derived a tensor basis for the photon/gluon self-energy and 
calculated the corresponding structure functions for an anisotropic system in which the distribution 
function is homogeneous and obtained from an isotropic distribution function by the rescaling of 
only one direction in momentum space 

/(p, x) = N{0 /iso (pv'l + e(v-n)2^ , (2) 

where N{(^) is a normalization constant, v is the particle 3-velocity, n specifies the anisotropy 
direction, and — 1 < ^ < cxd. Note that the distribution function given by Eq. is obtained in 
the limit ^ ^ oo assuming that Eq. ^ is normalized in the same way as Eq. 

Using classical kinetic field theory we were then able to numerically determine the photon /gluon 
self-energy structure functions for this tensor basis in the entire complex energy plane for arbitrary 
^. In this paper we will use these structure functions to determine the leading-order (collisional) 
energy loss of a heavy fermion propagating through an anisotropic QED plasma. This calculation 
will allow us to determine the dependence of the collisional energy loss on the angle of propagation 
On, the parton velocity v, the coupling constant e, and the temperature T. We will show that 
for large anisotropies and velocities the directional dependence of the heavy fermion energy loss is 
large and could therefore lead to a significant experimental effect. During the development we will 
also discuss some technical details related to the cutoff dependence of the isotropic and anisotropic 
results. 

The organization of the paper is as follows. In Sec^Jwe review some of the notation and results 
from Ref. In Sec. IlII Al we calculate the contribution to the energy loss coming from the 

exchange of soft photons with momenta on the order of eT. In Sec. IIIIBl we discuss complications 
which could arise due to the presence of unstable soft photonic modes and show that the soft 
energy loss calculation is safe. In Sec. lIIlCl we calculate the contribution to the energy loss coming 
from the exchange of hard photons with momenta on the order of T. In Sec. IIVI we combine the 
soft and hard contributions to obtain the final isotropic and anisotropic results. In Sec. [V] we list 
the limitations of the treatment presented here and provide a short description of how the result 
contained here can be extended to QCD. 



II. SETUP 



In this section we will review some of the findings from our previous paper which we will 
use in the subsequent sections. 



A. Tensor Basis 

For anisotropic systems with only one preferred direction we construct a basis for symmetric 
3-tensors that depends on a fixed anisotropy 3-vector n* with = 1. We first define the projection 
operator 

j^ij ^ §ij _ k^k^/k^, (3) 
and use it to construct n* = A'^n^ which obeys h ■ k = 0. With this we can construct the tensors 

B'^ = k'y/k"^ (4) 
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D'^ = k'n^ + k^n\ (6) 

Any symmetric 3-tensor T can now be decomposed into the basis spanned by the fom' tensors 
A,B,C, and D 

T = aA + 6B + cC + dD . (7) 
B. Self-energy structure functions 

Converting the result of Ref. from QCD to QED the spacehke components of the high- 
temperature photon self-energy for particles with soft momentum {k ~ eT) can be written as 

where the electron distribution function /e(p) is, in principle, completely arbitrary. In what follows 
we will assume that /e(p) can be obtained from an isotropic distribution function by the rescaling of 
only one direction in momentum space. In practice this means that, given any isotropic distribution 
function /e,iso(p), we can construct an anisotropic version by changing the argument of the isotropic 
distribution function 



/e(p) = N{i) /e,iso (^Vp' + ^(P • j ' (9) 

where N{£^) is a normalization constant, n is the direction of the anisotropy and ^ > — 1 is a 
adjustable anisotropy parameter. Note that ^ > corresponds to a contraction of the distribution 
along the n direction whereas — 1 < ^ < corresponds to a stretching of the distribution along the 
n direction. 

The normalization constant, N[£^), can be determined by normalizing the distribution function 
to a fixed number density for all values of ^ 

^ fe,Uv) = ^ /e(p) = N{£,) ^ /e,iso (yWOP^) ■ (10) 

By performing a change of variables to p 

p2=p2(l+C(vn)2) , (11) 
on the right hand side we see that the normalization condition given in requires that 



A^(0 = vr+e- (12) 

Making the same change of variables (|11|) also in Q it is possible to integrate out the \p\- 
dependence giving 



where 

^ - ^ dpp -r- 

Tx^ Jo dp 



^1 = -ZtL dpp'^^^^^- (14) 



Note that in the analysis which follows we will assume that /e,iso(p) = it-f{p) = (exp(p/T) + 1) ^, 
in which case mj^ = 
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We can then decompose the self-energy into four structure functions 

n = aA + /?B + 7C + JD , 
which are determined by taking the following contractions: 

k-n-k = (3, 



(15) 



n n k 
n • n • n 

Trn 



n^(Q + 7) , 
2a + /3 + 7 . 



(16) 



All four structure functions depend on mo, oj, k, ^, and k • n. In the limit ^ ^ the structure 
functions a and P reduce to the isotropic hard-thermal-loop self-energies and 7 and 6 vanish 



lim a(K) 



2 



lhn/3(i^) = _ni(K) + O(0, 



lim ^(K) 
lim 5(K) 

4^0 



O{0 

o{0 



with 



Ut{K) 



^ 1 
777 log 

2k UJ 



2ujk 
UJ + k 



A:^ u + k 



k 



1 



(17) 

(18) 
(19) 



The C'(^) terms in l|17|) were determined analytically in Ref. ^ 

With these structure functions in hand we can construct the propagator A*-'(/C) using the 
expressions from the previous section 



A(i^) = Aa [A - C] + Ag [{k^ - ^2 ^ a + 7)B + (/3 ■ 



■ UJ 



2^C 



where 



k^ — uj" + a 
u2 , ,2 



uj'^ + a + -f){l3 -uj'') -k^n^d"^ 



(20) 



(21) 
(22) 



III. ENERGY LOSS CALCULATION 

In this section we compute the collisional energy loss of a heavy fermion propagating through an 
electromagnetic plasma which has an anisotropic momentum space electron distribution function. 
The starting point for the calculation is the same as in the isotropic case In these papers 

Braaten and Thoma calculated the energy loss of a heavy fermion for both QED and QCD. Here 



Note that the normalization used in Ref. ll]J| is different than the normalization used here. The weak- 
anisotr opy^ expansions for the normalization contained here can be obtained from the expressions contained 
in Ref. lOj by taking m|j ^\ + ^ m|) and re-expanding in ^. 
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we will concentrate on QED but the extension to QCD is relatively st raig htforward. Note that 
there were previous calculations of the in-medium partonic energy loss jld. [l5|: however, the 
calculation of Braaten and Thoma was the first to present a systematic method for performing the 
complete order g calculation of this quantity. 

The technique used by Braaten and Thoma was to consider independently the contributions 
from soft photon exchange {q ^ mo ^ eT) and hard photon exchange {q ~ T). Each of these 
quantities is divergent with the soft contribution having a logarithmic UV divergence and the hard 
contribution having a logarithmic IR divergence. In their paper Braaten and Thoma found that 
if an arbitrary momentum scale q* was introduced to separate the hard and soft regions then the 
dependence on this arbitrary scale vanished when the contributions are combined and the result 
obtained is finite. 

Below we will use the same technique as Braaten and Thoma but we will show that if the 
resulting integrals are treated more carefully then the final result obtained by adding the soft and 
hard contributions does depend on q* . However, in the weak-coupling limit e <C 1 the dependence 
on q* is small as long as m/j ^ g* ^ T and the Braaten-Thoma result is reproduced. Since the 
full result does depend on q* we will need a prescription for fixing it. Here we will fix q* using the 
"principle of minimal sensitivity" which in practice means that we minimize the energy loss with 
respect to q* and evaluate at this point. We can then obtain a measurement of our theoretical 
uncertainty by varying q* by a fixed amount around this point. In the next three subsections we 
derive integral expressions for the soft and hard contributions to the energy loss and present some 
results for these in the limit of small anisotropies. We also discuss complications which could arise 
due to the presence of plasma instabilities. 



A. Soft Part 

For calculating the soft energy loss one can use classical field theory methods. The classical 
expression for parton energy loss per unit of time is 

(^) ft " / ' ^'"'^^^^ ' ^^^^ 

where X = (t,x), and Jgxt is the current induced by a test fermion propagating with velocity v: 

Jext((3) = (27r)ev(5(a; - q • v) , (24) 

with Q = (w, q). Using 

^Ld(Q) = (a'^(Q) - ^i{Q)) JiAQ) , (25) 

where Aq-' is the free propagator, and Fourier transforming to X we obtain 

EUX) = ie| ^ (q . v)(A^^(Q) - A^^(Q))^^e^('i-('i-)*) . (26) 
The above equation allows us to write Eq. H23|) as 
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The propagator can be expanded in our tensor basis and the contractions on the right-hand-side 
give 

v'B'^v^ = (q-vf, 
v'&^v^ = (n•v)Vn^ 

vW^^v^ = 2(q- v)(n- v) , (28) 
so that together with dW/dx = v~^dW/dt we obtain 



da; /soft V 



1 



cp- — ufi 



(n • v)^ 



^ 2 ~ ~2 



(27r)3 

2 l~ \2 

^{q^ - a;2 + a + 7) + (/? - u;^)!^?^ _ 25u{iv ■ v) 



q^ 



+ 



uj{q'^ — uj^) 



i4(,2_,2)_,2(n^ 



(29) 



with a; = q • V. Performing some algebraic transformations and scahng out the momentum gives 



V"d^Aoft " V (2^^(T^a;2) 



—a 



(g2 — q^u"^ + a) 

q^A + B 



I 2 -2 
[V — U 



(n • v)^ 



q'^C + q^V + S 



(30) 



where 



C 
V 
£ 



V? 



-(a + 7) - 2(j(l - w2)(n- v)(5 , 



B = {{a + -i)l3-iiH'^) (1-w 



.2 (a-v)^ 

~9 / ' 



-a;2(l-(2;2)^ 
{a + 7)/3 - n^^^ , 



(31) 



with a) = q • V and 5 = q6. With all the momentum dependence made explicit we can now perform 
the q integration to obtain 



/dW\ 



e" f dQ.„ u) 
— Im ' 



soft V 



-a — rrrr^ In ■ 



where 



(27r)3 (1 - (£;2) " 2(1 - w^) a 

+F{q*) - F(0) 



AV - 2BC ^ ^JV^ - ACS + V + 2Cq' 



-An 



4C£ VV^ - 4C£ -V- 2Cg2 
and we introduced a UV momentum cutoff q* on the q integration. 



(32) 



(33) 
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10"' 10"' 1(f 10' irf 

q* /mo 

FIG. 1: Isotropic energy loss of the soft part 
(full line) function of q* /rriD for v = 

0.5 compared to the result from Braaten and 
Thoma (dashed line). 
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FIG. 2: The contribution — {dW / dx) ^^^^ as a 
function of q* /m£, for v = 0.5. 



1. Small-^ Limit 



Using a linear expansion in the limit of small ^ one can obtain analytic expressions for the 
structure functions a, (3, 7, and 5 as given in Ref. A subsequent expansion of Eq. to 

linear order in ^ gives 



V dx / soft, small- 



/dT^\ 



soft,iso 



/dVF\ 



soft,i 



(v- n)2 /dVF\ 
V dx / 



soft, ^2 



(34) 



"where 



dx 



Im 



soft,iso 



dCj 

(2^ 



■ UJ 



2(1 



2^2 



Ht In 



(1 - (2;^)^*^ + nr 



2 



*2 



Hi 



-Hi 



(35) 

is the isotropic result, "which - assuming (f ^ m/) - corresponds to the result from Braaten and 
Thoma "when inserting the explicit form of W-Tifii) and Jlj^iCS) from Eqs. p8j) and (|19|) . The form 
of the contributions (dVF/dx)gQf^ and (dTy/dx)gQ£^. resembles that of (dTy/dx)gQ££ jg^, but since 
they consist of many more terms than the isotropic contribution "we refrain from giving them here 
explicitly. 

Note that by scaling out the Debye mass niD from the self-energy functions, the contribu- 
tions (dTy/dx)g(,f^ jgQ, (dTy/dx)gQ£^ and (dVF/dx)gQ£^ depend on the particle velocity the ratio 
g*/m£), and an overall multiplicative factor of e^m^j only. The direction of the heavy fermion 
enters the whole result Eq. H34() only through the explicit term (v • n)^/f^ in the small ^ limit. 
Note that, in principle, the presence of unstable modes "would make the g-integration divergent; 
ho"wever, this is not the case for the soft energy loss as "we will discuss in Section [ill Bl 



2. Behavior of the Soft Part 

The isotropic result for the soft part Eq. (|35|) (corresponding to the one obtained in Ref. [l^ ) 
is shown in Figure Q together with the Braaten-Thoma result from Ref. ^3] for v = 0.5. As can be 
seen, the results are identical for large q*/mD, while for small q* /rriD the Braaten-Thoma result 
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becomes negative whereas the full result obtained by numerical integration of Eq. I|35|) is positive 
for all values of q* /niD- 

In Figured we show the function — [dW / dx) ^^^^ which controls the directional dependence 
of the soft part of the energy loss for small ^. As can be seen from this Figure the function is 
negative for small q* /m£, but becomes positive at a velocity-dependent value oi q* /m£,. The value 
of q*/mD where {dW / dx) ^^^^ = is of special interest here since it signals the point at which 
the directional dependence of the energy loss changes. The value of q* /rriD obtained using the 
principle of minimal sensitivity is a function of the velocity and coupling constant and therefore we 
expect a non-trivial dependence on these parameters. For ^ > this means that as the coupling 
constant is increased from zero for fixed v that at first the energy loss will be peaked along n but 
will eventually become peaked transverse to n. The value of the coupling at which this change 
in the directional dependence occurs increases as the velocity of the fermion increases. For the 
exact expression given by Eq. (|32() the results obtained in the small ^ limit hold qualitatively but 
quantitatively the predictions differ considerably once ^ becomes large. 



3. Large-^ Limit 

In the limit ^ ^ oo one finds that the anisotropic distribution function turns into a specific case 
of Eq.dU), 

^lim ^l + inp {p^jl + C(P • n)2^ = 5(p • n) j dx np (pVTTx^) , (36) 

where the remaining integral representin g F (p \) in Eq.(^ can be evaluated analytically to be a 
sum over Bessel functions. Following Ref. |11| it is then possible to evaluate the structure functions 
analytically. The result for the soft part of the energy loss is then given by Ea. H32|) with the general 
structure functions replaced by their large-^ expressions. For large values of ^, the general result 
for the soft part converges towards the limiting result, as it should. 



B. Dynamical Shielding of Plasma Instabilities 

As shown previously 0, 0, 0, 0, 0, 11 1 in systems in which the distribution function is 
anisotropic in momentum space there exist unstable modes which could potentially make physical 
quantities incalculable within a perturbative framework However, we will show in this sec- 
tion that in the case of the energy loss the singularities induced by the presence of these plasma 
instabilities are "shielded" and are therefore rendered safe. To see this one need only consider the 
contribution from the propagator to Eq. H3U|) which is schematically given by 

~ ^™ / / qdq ^ 2^2-^ — V • 

Vdx/A.soft J J (1 - w^) (g^ - g^w^ + a) 

The possibility of singularities arise because in the static limit the structure function a is negative- 
valued which can result in singularities along the integration path which are unregulated. However, 
in the limit uj ^ the structure function a has the form 

\ima{uj,q) = M'^{-l + iDi:j) + 0{uj'^) , (38) 

where M and D depend on the angle of propagation with respect to the anisotropy vector and the 
strength of the anisotropy. 
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0.6 
^ 0.4 



' 10"^ 10"' 1(f lrf 1(f 

FIG. 3: Coefficient D as a function of ^ for FIG. 4: Coefficient as a function of 0„ for 
9n = 1.5. i = 100. 

Ignoring the 0{lo) contribution above we see that a singularity in the integrand can arise at the 
point q = M when — > 0. However, including the 0{uj) contribution we see that the combination 
of the power of cD in the numerator and the presence of a contribution proportional to c!) in a 
together render the singularity finite at this point as long as D is non-vanishing 

hm hm — — — — > — — — > — . (39) 

d;^o g^M - q^(:b^ + M^i-l + iDib) M^{Qj + iD) M^D ^ ' 

Therefore, the g-integration can safely be performed. Of course, the integrand may change rapidly 
in this region so it should be treated carefully when numerically evaluating the integral. Note that 
this is similar to dynamical screening of the magnetic sector of finite temperature QCD. In fact, 
in the isotropic limit, the coefficient D approaches the well-known value of —mj^ which is due to 
Landau damping. Unfortunately for anisotropic systems it is difficult to prove that the coefficient 
D is non-vanishing for all values of the propagation angle and anisotropy strength. However, we 
can evaluate D analytically in the weak-anisotropy limit and numerically for general ^ and in both 
cases we find that D is non-vanishing. 

In the weak-anisotropy limit we can derive analytic expressions for all of the structure functions 
as shown in Ref. To linear order in ^ the structure function a becomes 



lim lim a = mn 



-|(1 + cos26'„) - f 1 -h ^(1 + cos26l„) 



where 0„ is the angle of propagation with respect to the anisotropy vector n. Therefore, in the 
weak-anisotropy limit D is non-vanishing for all ^ > 0. For ^ < 0, the mode is stable so there 
is no singularity to be concerned about. At higher orders in the weak-anisotropy expansion the 
picture becomes slightly more complicated but is still consistent with D being non-vanishing. 

For large values of ^ we can not rely on a weak-anisotropy expansion so we have to resort to 
numerical determination of D. In Fig. El we plot the dependence of D on ^ for On = 1.5 and in 
Fig. |1] we plot the dependence of D on for ^ — 100. As can be seen from these Figures, D is 
non- vanishing for all values of ^ and 9n shown. Beyond what we have plotted, we have calculated D 
for numerous values of ^ and 9n and in every case we find that D is non- vanishing. We are therefore 
reasonably confident that D is in general a non-vanishing quantity and therefore the singularities 
which could come from the unstable modes are dynamically shielded and thereby rendered safe in 
the energy loss calculation. Note that similar arguments can be made to show that the singularities 
coming from the Aq contribution to Eq. 1)30^ are also dynamically shielded. 
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FIG. 5: Tree-level Feynman diagrams for the scattering processes e~ fj, — > e~/u and 7/x jfi. 
Note that the second two diagrams cancel with each other so that in QED only the first diagram 
contributes to the hard energy loss. 

C. Hard Part 



The expression for the hard contribution to the energy loss for an arbitrary electron distribution 
function /e(p) is obtained by summing the tree level diagrams shown if Figure El Assuming 
V ^ T/E, performing the Dirac traces, and summing over spins the result can be reduced to 



\ dx / hard 



Aire^ f d^k /e(k) f d^k' 1 - /e(k' 



(2^)3 k J (27r) 



-6{u; — V • q)0(g — q*) 



,2X1 



2(A:-v-k)(A;'- vk') + 



(41) 



where uj = k' — k and q = k' — k, and we have introduced an infrared cutoff (f on the q integration. 
The term involving the product /e(k)/e(k') vanishes since the integrand is odd under the inter- 
change of k and k'. Redefining the origin of the k' integration so that it becomes an integration 
over q we have 

"4 " d^k /e(k) 



/ hard 27r2w } i^Tx)^ k 



q dq dO. 



5{u) — V • q) 

|q + k| 



UJ 



2(A: - V • k)(cj A; - V • k - V • q) 



, (42) 



where now iv = \q + k\ — k. Choosing v to be the z-axis for the q and k integration we can rewrite 
the delta-function as 



(5(|q-hk| -k-vq) 



,5((/)q-(/)o)G(A; + vq)2|q + k| 



qsJAk'^ s\T? 6k sin^ 6q - {q{l - cos^ Bq) + 2cos9g{kcos9k - kv)Y 
where (/>o is the solution of the equation 

q{l — cos^ 0g) -|- 2 cos Qq(k cos — ^v) 



cos((/)o - (\>k) 



Ik sin Qk sin Qn 



(43) 



(44) 



The <^q integration is then straightforward and we obtain 



da; / hard 27r2t; j (27r) 



S'k /e(k) 

g ~~I — / Qdq d cos dq 4Q{k + vq cos 9q 



e(4A;2 sin^ 9k sin^ Oq - {q{l - cos^ 9q) + 2 cos 9q{k cos 9k - kv)f) 
4A;2 sin^ 9k sit? 9q - (g(l - v"^ cos^ 9q) + 2 cos 9q{k cos 9k — kv)Y 

2(/t-vk)2 + 1^(0,2 _g2) 



(a;2 



(45) 
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where uj = vqcosOq and we have included a factor of 2 because of the symmetry (/>o <-> 27r — 
Setting fe,\so{k) = np^lt) and scahng k = qz one can perform the q integration to obtain 



dW 



dx 



hard 



{q*fT^VTT^ 
{2tt)^v 



zdz J d cos Ok 



d(t)kFi{q*z^Jl + ^{n^ sinOk cos (pk + n^ cos Ok)'^) 

9(4^2 sin^ Ok sin^ Og - {I - cos^ 6*5 + 2 cos 6lgz(cos Ok - v)f) 



fi cos 9q Q{z + V COS 0q) 

-1 ' Uz^ sin2 sin^ 0q - {1 - cos2 0^ + 2 cos 9qz{cos Ok - v))^ 



V COS 9„ 



{v"^ cos^ 6*^ — 1) 
where q* = q* /T, 



2^2(1 cos 6*^)2 + 



1 - 



(W^ COS^ Oq-1) 



Fi{x) 



X In (1 + exp (— x)) — Li2(— exp (— x)) 



(46) 



(47) 



and TT-z = n • V with 1 = + n^. To obtain the final result for the hard contribution to the energy 
loss the remaining integrations have to be performed numerically. 



1. Small-^ Limit 

For small ^ we use Eq. H42|) and scale out the ^-dependence using Eq. ® and substituting 
^ /^^(l + .^(k • n)2) and q^ ^ (/^(l + ^(k • n)^). Performing a linear expansion in ^ then gives 



fdW\ _ 

V /hard,small-C '^T^'^V Jo 



where 



/e,iso(^) 



q* 



9i{q) - C 92{q) 



51(g) 



(48) 



9iiq) 



92{q) 



dn. 



dn. 



dVlk 2 '^('^ ~ ■ ^) 



UJ 



4tt 



|q + k| 
■q'{k.n) 



2(A;- vk)2 + 



1 



d^k 2/r \2^(.^ - V • q) 



UJ 



2(A:- vk)2 + 



1 



■ (49) 



An |q + k| uj^ — q^ 

Taking q as the z-direction makes the d^lk integration straightforward but still algebraically 
intensive for the 32 contribution because of the extra dependence on k • n. Next we rotate the 
coordinate system so that v becomes the new z-direction and n lies in the x — z plane for the d^lq 
integration; the polar integration is then also straightforward whereas for the azimuthal integration 
the limits coming from the delta function in Eq. (|49|) are somewhat non-trivial. Nevertheless, by 
scaling k = qz one can perform the q integration to obtain 



dW \ 

/ hard, small— ^ 

where q* = q* /T, 



e^{q 



*\2 



27r2^ 



dz 
2^ 



Fiiq*z) 



gi{l) = — dujQ{v 

ZV J\l-z\-z 



W^) U> 



5i(i)-e 52(1) 



— + 3z(z + uj) 
4 



9i{l) 



^fiso{q*z)g2{l] 



(50) 



(1-.2) 



z{z + Uj) 



1 



(51) 
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FIG. 6: Isotropic energy loss of the hard part 
(full Hne) function of q* for t; = 0.5 com- 
pared to the result from Braaten and Thoma 
(dashed line). 



10-^ 



FIG. 7: The contribution - 
function of q* for v = 0.5. 



{dW/dx] 



hard,52 



as a 



and the exact form of §2 is similar to gi but involves many more terms so we have refrained from 
giving an explicit expression. The integration over uj can then be performed analytically using 

r{l+v)/2 



dz 



l—z\—z 



dz 



duj + 



1-22 



dz 



{l+v)/2 



du) . 



(52) 



l{l-v)/2 

but the remaining integral over z has to be performed numerically. The final result to linear order 
in ^ can then be written as 

(v-n)2 fdW 



fdW\ 



hard, small— 5 



/dW\ 



hard,iso 



/dW\ 



hardj^l 



dx 



hard,52 ^ 



(53) 



where idW/dx)^^^.^-^^^ is given by Eq. (jSU)) with ^ = 0. 



2. Behavior of the Hard Part 



In Figure ini we show the result for the isotropic (^ = 0) result for the hard part compared to the 
Braaten-Thoma result from Ref. ;12]. Similar to what we found for the soft part, the two results 
are identical for very small q*, while for large q* the Braaten-Thoma result becomes negative and 
our result is positive for all values of q*. 

In Fi gure [3 we plot the function — (^dW/dx)^^^.^ a which controls the directional dependence 
of the hard part of the energy loss. As was the case with the soft part, idW/dx)^^^^^ = for a 
velocity-dependent value of q*. Again, in practice, this means that the directional dependence of 
the hard energy loss changes as the coupling constant is increased. 



3. Large-S, Limit 



Using the explicit form 1)36(1 for /g one can do the dcj))^. integration in Ea. (|45() by rewriting 
r/f M 5{(j)k-4>o) n cose^fcn^ 

(5(kTi) = — cos0o = — —f, • (54) 

Vn^ - cos^ 6*^ sm BkU^ 
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FIG. 8: Energy loss as a function of q* for v = 
0.5 and e = 0.01. The dashed hne corresponds 
to the result from Braaten and Thoma. 
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FIG. 9: Energy loss as a function of q* for v = 
0.5 and e = 0.5. The dashed line corresponds 
to the result from Braaten and Thoma. 



After scaling cos 6k cos 9k and doing the g-integration the remaining integrals take the form 

'dW\ ?>2e\q*fT'^ ^ f^^^,^, n ^, ^ 1 ^ n 

— — — ^ / zdz / dxFi[q z\/l + x^) / dcosdk—r= / dcosO„ 

{2tt)^v Jo Jo J^i VI - cos^ 6*^ J-i 

6(42^(1 - nl cos^ 9k) sin^ 9g - {I - v"^ cos^ 9q + 2cos.9qz{nx cos9k - v))'^) 



) hard 

0(z + V COS 9qj 



V COS I 



(f 2 COS^ 9q — 1)2 



4^2(1 — cos^ 9k) siv? 9q — (1 — v'^ cos^ 9q + 2 cos 9qz{nx cos ^a: — v))'^ 
2z^(l-t;n^cos6lfc)^+ (i;^ cos^ (9g - 1) \, (55) 



which after evaluation then give the result for the hard part of the energy loss in the large-^ limit. 



IV. NUMERICAL EVALUATION 



A. Isotropic Case 

In the isotropic case, the total collisional energy loss is obtained by adding Eq. and ij^O]) 
with C = 0, 

(dW\ (dW\ (dW\ 



+ 



V dx ; 



(56) 



*l2^ / iso.soft V / isOjhard 

In general the isotropic energy loss is a function of the electromagnetic coupling e, the velocity of 
the particle v, the temperature T, and the momentum separation scale q* . 

In principle q* should be an arbitrary quantity in the range e <C g*/r <C 1 which is possible in 
the weak-coupling limit. Choosing e = 0.01 and plotting —dW/dx as a function of q* = q* /T at 
fixed V and T we find the result shown in Figure El From this Figure we can see that the energy 
loss develops a plateau where the value is essentially identical to the result obtained by Braaten 
and Thoma; however, outside the plateau the function rises logarithmically with q* . For higher 
coupling (e ^ 0.5) the plateau shrinks rapidly to a minimum which can be seen in Figure 1^1 We 
fix q* by minimizing the energy loss with respect to q*. 

In Figure [TUl we compare our result with the result obtained by Braaten and Thoma for e = 0.3. 
The Braaten-Thoma result is shown as a dashed line and our result is shown as a grey band. The 
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FIG. 10: Isotropic energy loss as a function of the electron velocity v for e = 0.3. Dashed line is the 
Braaten-Thoma result. The grey band indicates the variation of the result obtained by varying q* 
by a factor of two around the central value. Inset is an enlarged version of the region 0.3 < f < 0.7. 



band corresponds to the variation we obtain in our final result when varying the scale q* by a 
factor of two around the central value (minimum). Note that the choice of varying by a factor of 
two is completely arbitrary. We could have easily chosen a smaller variation, however, we have 
chosen to be conservative here and vary q* by a factor of two. As can be seen from this Figure, for 
e = 0.3, the corrections to the Braaten-Thoma result are small; however, for larger couplings the 
band obtained by varying q* by a factor of two can be quite large. The size of these bands gives 
us an estimate of the theoretical uncertainty present in the calculation. 



B. Anisotropic Case 

In the general anisotropic case we have to numerically evaluate the integrals contained in 
Eqs. (|32p and (|46|) and add these to obtain the total energy loss 

V dx / V dx / soft V dx / hard ' 

In addition to v and e the anisotropic energy loss also depends on the anisotropy parameter ^ and 
the angle of propagation with respect to the anisotropy vector On- Note that for every value of e, 
^, and On there is a non-trivial consistency check between the hard and soft contributions, namely 
that the coefficients of log q* have to be equal and opposite for small-^* and large-g* , respectively. 
In every case that we have examined we find that these coefficients match to a part in 10~^ which 
can be further improved by increasing the target integration accuracy. 

As in the isotropic case we find that for small coupling there is a plateau for e <C <?* ^ 1 and 
our final result is obtained by finding the minimum of the energy loss as a function of q* . This 
corresponds to the point at which the result is minimally sensitive to q* . As before, we can then 
obtain an estimate of the uncertainty in our prediction by varying q* by a factor of two around the 
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FIG. 11: Scaled energy loss function of 
the angle of propagation with respect to n for 
V = 0.3 and = 1. 



^1^ 



/ 

y □ e = 0.05 

LJ e = 0.3 

n e = 1.0 



FIG. 12: Scaled energy loss as a function of 
the angle of propagation with respect to n for 
V = 0.3 and ^ = 10. 




FIG. 13: Scaled energy loss as a function of 
the angle of propagation with respect to n for 
V = 0.5 and ^ = 1. 
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FIG. 14: Scaled energy loss as a function of 
the angle of propagation with respect to n for 
V = 0.5 and ^ = 10. 



point of minimal sensitivity. Note that using this method the central value of q* always serves as 
a lower-bound on the prediction. 

In Figures I11I16I we plot the ^^-dependence of the resulting anisotropic energy loss for v = 
{0.3,0.5,0.7}, e = {0.05,0.3,1.0}, and ^ = {1,10} scaled by the corresponding isotropic energy 
loss. For e = 0.05 we see that for all three values of the velocity the energy loss is larger along the 
direction of the anisotropy then transverse to it. For the physical value of electromagnetic coupling 
constant, e = 0.3, we find that directional dependence of the energy loss depends on the velocity 
and the strength of the anisotropy. For v = 0.3 and ^ = 1 the result is almost independent of the 
direction of propagation but if the anisotropy is stronger = 10) then the energy loss is less in 
the direction of the anisotropy than transverse to it. Note that this is consistent with the results 
obtained in the small-^ limit in Sections IIII A 21 and IIII (T2l For larger velocities we find that the 
energy loss along the direction of the anisotropy is higher than transverse to it; however, if the 
coupling constant is increased to a large enough value then the trend will reverse as was the case 
for smaller velocities. Note that for v ^ 0.5 the coupling constant required to reverse the trend is 
extremely large. 
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FIG. 15: Scaled energy loss function of 
the angle of propagation with respect to n for 
V = 0.7 and ^ = 1. 



FIG. 16: Scaled energy loss as a function of 
the angle of propagation with respect to n for 
V = 0.7 and ^ = 10. 



In the most extreme case shown here in Fieure lT^ for which v = 0.7 and ^ = 10 we see that for all 
values of the coupling constant shown the energy loss varies by 20-25% depending on whether the 
fermion is propagating along the direction of the anisotropy or transverse to it. For larger values 
of the velocity and anisotropy parameter the directional dependence of the energy loss increases. 

In Figure El we plot the dependence of the energy loss normalized to the isotropic result for 
^ = 1, e = 0.3, and 6'„ = {0, 1.5} as a function of the velocity v. From this Figure we see that 
for fixed e and ^ the angular dependence (forward-peaked versus transversely-peaked) changes as 
the velocity is increased with the energy loss being transversely-peaked for small velocities and 
forward-peaked for large velocities. For the largest velocity shown in this Figure we see that for 
^ = 1 the effect is on the order of 10 %. 

In Figure ITHl we plot the dependence of the energy loss for v = 0.5, e = 0.3, and 9n = {0, 1.5} 
as a function of ^. Additionally, we have included results obtained by taking the large-^ limit 
by lines on the right hand side of the plot frame. As can be seen from this Figure for ^ > we 
find that the difference between the forward and transverse energy loss saturates at approximately 
25% in the large ^ limit. For ^ < we see that for both angles shown the energy loss approaches 
zero but the difference between the forward and transverse energy loss is still quite large at = 
—0.99. Therefore, it is possible that the presence of momentum-space anisotropics in the electron 
distribution function could lead to a rather significant experimental effect. 



V. CONCLUSIONS 

In this paper we have derived integral expressions for the collisional energy loss of a heavy 
fermion propagating through a QED plasma for which the electron distribution function is 
anisotropic in momentum space. We then numerically evaluated the resulting integrals and studied 
the dependence of the heavy fermion energy loss on the angle of propagation, degree of momentum- 
space anisotropy, and coupling constant. We have shown that the techniques used in the isotropic 
case can be straightforwardly extended to the anisotropic case. We have also discussed how prob- 
lems due to unstable soft photonic modes are avoided in the case of the energy loss. 

As a side result we demonstrated that in the isotropic case the canonical result from Braaten and 
Thoma has a correction which comes from the dependence of the result on the scale introduced 
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V 1+e 
FIG. 17: Energy loss for angles 6n = {0, 1.5} FIG. 18: Energy loss for angles 9n = {0, 1.5} 
as a function of v for ^ = 1 and e = 0.3. as a function of for v = 0.5 and e = 0.3. Lines 

on right hand side of the plot frame indicate 
results obtained in the large-^ limit. 



to separate the soft and hard contributions to the energy loss. The residual dependence on the 
separation scale was then used to estimate the theoretical uncertainty of the resulting calculation 
for both the isotropic and anisotropic case. 

Our final results indicate that for anisotropic QED plasmas there can be a significant directional 
dependence of the energy loss for highly-relativistic fermions if there is a strong momentum-space 
anisotropy present in the electron distribution function. It should be mentioned that we have 
assumed here that the heavy fermion is infinitely massive. In the realistic case when the fermion 
mass is on the order of the muon mass the result obtained here should be reliable for velocities 
0.6 < I- < 0.95 [3. 

We have concentrated here on the collisional energy loss in QED. In order to extend the result 
to QCD there will be a modification of the hard energy loss due to the fact that the fermion- 
boson scattering diagrams do not cancel against one another in this case. Besides this, the only 
modification required will be to adjust the Debye mass so that it corresponds to the QCD Debye 
mass. One complication will come from the fact that the quark and gluon distribution functions can, 
in general, have different momentum-space anisotropics; however, there should be an independent 
matching of logarithms coming from the quark and gluonic sectors. This work is currently in 
progress. 



Acknowledgments 

M.S. and P.R. would like to thank A. Ipp, A. Rebhan, and M. Waibel for discussions. M.S. was 
supported by the Austrian Science Fund Project No. M689. P.R. was supported by the Austrian 
Science Fund Project No. P14632. 



[1] K. Adcox et al. (PHENIX Collaboration), Phys. Rev. Lett. 88, 192302 (2002). 
[2] J.L. Nagle et al. (PHENIX Collaboration), arXiv: nucl-ex/0209015. 



18 



[3] L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rcpt. 100, 1 (1983); A. H. Mueller and 
J. w. Qiu, Nucl. Phys. B 268, 427 (1986); J. P. Blaizot and A. H. Mueller, Nucl. Phys. B 
289, 847 (1987); L. McLerran and R. Venugopalan, Phys. Rev. D 49, 223 (1994); 49, 3352 
(1994); 50, 2225 (1994); 1. Balitskii, Nucl. Phys. B 463, 99 (1996); J. Jalilian-Marian, 
A. Kovner, L. McLerran and H. Weigert, Phys. Rev. D 55, 5414 (1997); Yu.V. Kovchegov 
and A.H. Mueller, Nucl. Phys. B 529, 451 (1998); Yu.V. Kovchegov, Phys. Rev. D 60, 
034008 (1999); 61, 074018 (2000). 
[4] R. Baier, A.H. Mueller, D. Schiff, and D.T. Son, Phys. Lett. B 502, 51 (2001). 
[5] E.S. Weibel, Phys. Rev. Lett. 2, 83 (1959). 
[6] S. Mrowczyhski, Phys. Lett. B 314, 118 (1993). 
[7] S. Mrowczyhski, Phys. Rev. C 49, 2191 (1994). 
[8] S. Mrowczyhski, Phys. Lett. B 393, 26 (1997). 
[9] J. Randrup and S. Mrowczyhski, arXiv: nucl-th/0303021, (2003). 
[10] P. Romatschke and M. Strickland, arXiv: hep-ph/0304092, (2003). 
[11] P. Arnold, J. Lenaghan, and G.D. Moore, arXiv: hep-ph/0307325, (2003). 
[12] E. Braatcn and M. Thoma, Phys. Rev. D 44. 1298 (1991); 44, R2625 (1991). 
[13] J.D. Bjorken, Fermilab Report No. PUB-82/59-THY, (1992). 
[14] M.H. Thoma and M. Gyulassy, Nucl. Phys. B 351, (1991). 
[15] S. Mrowczyhski, Phys. Lett. B 269, 383 (1991). 

[16] G.D. Moore, QCD Kinetic Theory (Lecture), Quantum fields in and out of equilibrium. 
University of Bielefeld, (2003). 



19 



